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SEA-AIR  EXCHANGE  OF  ENERGY  AND  MOMENTS  UNDER 
WELL-DEVELOPED  SEA  CONDITIONS:  THEORY  AND  EXPERIMENT 

V.E.  ZAKHAROV,  P.l.  AND  E.  KUZNETSOV,  CO-PI 

Ocotober  1,  1994-Sept0mber  30,  1995 

Our  work  in  1995  is  going  very  succesfully.  We  received  following  essential 
results: 


l.We  achieved  a  serious  progress  in  numerical  simulation  of  surface  waves 
on  free  surface  of  3D  deep  fluid  in  presence  of  gravity  and  surface  tension. 
At  the  moment, using  a  kind  of  spectral  code, we  can  perform  simulation  of 
256  X  256  spectral  modes.  In  this  problem  we  use  expansion  of  the  Hamil¬ 
tonian  up  to  terms  of  fourth  order  in  canonical  variables.  The  developed 
technique  allows  to  solve  a  lot  of  different  physical  problems.  So  far  we  used 
it  for  solving  a  problem  of  a  big  theoretical  importance.  By  direct  simula¬ 
tion  we  comfirmed  existence  and  stability  of  the  weak-turbulent  Kolmogorov 
spectrum  for  capilliary  waves.  The  result  is  submitted  for  publication  to  the 
"Physical  Review  Letters”  (Pushkarev,  Zakharov). 


2.We  obtained  important  results  in  the  theory  of  2D  fluid  with  a  free 
surface.  Using  conformal  mapping  to  half  plane  directly  in  the  principle  of 
Hamilton,  we  found  a  way  of  decsription  of  the  surface  which  is  very  suitable 
for  numerical  simulation  by  spectral  code.  The  simulation  can  be  done  in 
exact  equations  without  any  assumpions  on  a  level  of  nonlinearity.  The 
developed  numerical  scheme  allows  to  describe  strongly  nonlinear  prosesses 
like  wave-breaking.  By  separating  of  scales  we  found  an  exactly  integrable 
model  in  the  the  theory  of  deep  fluid  with  free  surface.  The  model  can  be 
applied  to  analytical  decsription  of  the  wave  breaking  and  to  the  theory  of 
nonlinaer  phase  of  the  Raileigh-Taylor  instability.  The  results  were  reported 
on  the  International  Symposium  on  Applied  Mathematics  in  honor  of  70  year 
of  M.Kruscal(Boulder,  Colorado,August  1995)  (Dyachenko,  Zakharov). 


3. We  perfomed  a  numerical  calculation  of  the  Kolmogorov  constant  in 
the  theory  of  turbulence  of  ideal  fluid  in  one-loop  approximation  in  Clebsh 
variables.  The  article  is  submitted  to  the  "Physica  D”  (Balk,  Pushkarev, 
Zakharov). 
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Abstract 


The  second  method  of  calculation  of  Kolmogorov  constant  is  based 
on  the  calculation  of  the  analytical  expression  for  matrix  element  over 
eight  angles  in  /i-space.  Once  this  expression  for  averaged  matrix 
element  S{k^ki^k2^k3)  is  known,  the  problem  of  calculation  of  the 
constant  is  reduced  principally  to  the  numerical  calculation  of  the 
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1  Calculation  of  the  collision  integral  on  the 
scaling  spectrum 

Consider  four-wave  kinetic  equation 

~dt  ~  I  +  ^1  +  ^2  +  k^)6{u  4-  Wi  -  u;2  -  ^3) 

+  ^)dhdhdh  (1) 

"-k  '‘'ki  ^'ks 

We  are  introducing  the  spherical  system  of  coordinates  in  A'-space,  choos¬ 
ing  one  of  the  axises  directed  along  one  of  the  diagonal  of  fourangle  k,  k^ ,  k2,  k^ 
(see  Fig.l)  ; 


(k  cos  6,  k  sin  6  cos  6,  k  sin  0  sin  (f)) 

{ki  cos  Oi,  ki  sin  6i  cos  4>i.  ki  sin  6i  sin  0,),  i  =  1,3 

We  are  assuming  further  that  spectrum  and  dispersion  are  invariant 
with  respect  to  rotations  in  A'-space.  After  averaging  over  four  angles  6  and 
four  angles  ^  kinetic  equation  (  1)  takes  the  form: 


dnk  _ _ 

dt  ~  k‘i 


-I  dk  /"  *^(^5  ?  ^2?  ^3) 


dk  dki  dk2  dk^ 

dijo  duJi  duj2  doj^ 


'^k2  '^k^ 


'  1 

- [_ 

^n.k  Uk, 


- - — ^  ^(u;  +  ijJi  —  UJ2  —  u?3)du;ida.’2da?3 


riko 


Here  the  function  S{k^  ki,k2,  ks)  is  the  result  of  averaging  of  matrix  ele- 
ment  over  angles  $  and  (f)  in  A'- space 


S{k,  ^1,  ^2,  ^3)  ==<  ki  k2-k  ks)  >  {kkik2ksY  ^  = 

=  j  +  ^2  +  ks)  sin  0  sin  sin  62  sin  dsd9d6id62d0s  x 

d(f>d(l)id(f>2d(f>s{kkik2ks)'^~^ 
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Brackets  <>  mean  an  averaging  over  /?-space  angles. The  domain  of  in¬ 
tegration  here  is  [0, 2;r]  for  angles  0  and  [— |,  |]  for  angles  (j). 

The  frequency  and  scaling  spectrum  are  defined  as 


oj,,  =  C{u)k°' 


Uk  =  Rk-^  = 


(2) 


where  the  functions  C(i')  and  A{i/)  are  defined  by  (...)  and  (...)  cor¬ 
respondingly.  The  expressions  for  them  are,  however,  unimportant  for  our 
purposes,  because,  as  we  will  see,  they  drop  out  from  the  final  expression  for 
the  Kolmogorov  constant. 

We  will  mention  here  only  the  relations  following  from  our  definitions: 


ax  =  1/ 

R  =  C-^A{iy) 

On  scaling  spectrum  (  2 )  the  kinetic  equation  becomes 


dn-k  f  A\ 


6{u>  +  iOl  -  U)2  -  0J3)dl0iduJ2<kO3 

Function  S{k,ki,k2,k3)  is  a  homogeneous  one  with  respect  to  variables 


Lo: 


1  1  1  3  j  1  1  1 

5'((ea;)«,(eu;/:j)a,(eu;fc2)-,(ew/,3)°)  =  e'*  S{uj^ 
where  the  pow'er  of  homogeneouty  is 

2/3 +  Zd 


7  = 


-  1 


a 


Using  the  property  of  gomogeneouty  of  the  function  S  one  can  ”pull-out” 
the  coefficient  C  out  of  the  integration  sign: 
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duk 

dt 


L  L 

^k2  ^  y 


—  u}%^  —  6{u  +  ui  —  LO2  —  uj3)dujiduj2duz 


Here  the  integration  is  carried  out  over  infinite  domain  (see.Fig.2),  which 
we  will  break  into  four  domains  Making  conformal  mappings 

of  domains  II, III, IV  into  domain  I  (see  [1])  we  get 


6(^UJ  —  ^2  —  i03^diO-[diJ2diOj 


where 


?/  =  -7  +  3x  -  3 

Being  rewritten  in  a  form  (  3),  kinetic  equation  explicitly  exibits  four 
stationary  scaling  solutions  [1].  One  of  them  corresponds  to  y  =  1  and 
defines  the  final  energy  flux  P  from  long  to  short  scales.  Rewriting  integral 
(  3)  in  dimensionless  form  on  can  get  the  expression  for  a  collision  integral 
on  scaling  spectrum  (  2): 


where 


duk 

dt 


Hy) 


i{y)  =  jj  (666)""(i + er-  ef  -  m 

(1  +  —  ^2  ~  d"  ~  ^2  “  ^z)dw\dl02dujz 

Here  I,-  =  (^)“  are  dimensionless  frequencies,  2  =  1,3. 
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Domain  of  integration  0  is  defined  from  the  condition  that  four  vectors 
k,ki,k2.,k3  can  form  a  fourangle  on  manifold 


(jJ  UJ I  —  u^2  4"  ^3 

In  other  words  domain  0  is  defined  from  the  conditions 

ki  k2  ks  >  k 
CJ  -|-  =  LO2  “t"  ^3 

Uk  =  C{v)k°‘ 

.  a;2  <  a;,  Wa  <  a; 

Rewriting  last  conditions  in  dimensionless  units,  we  get; 

^  {  (^2  +  ^3  —  1)  “  +  ^2“  +  ^3  >  1)  ^2  <  1?  <^3  <  1  } 

This  domain  of  integration  D  is  shown  on  Fig.3. 


2  Energy  flux  and  Kolmogorov  constant 


The  energy  flux  P  from  the  pumping  region  of  long  scales  to  the  damping 
region  of  short  scales  at  intermediate  ’’transparency”  region  of  frequencies 
o-’jt  is  defined  by  the  equation 


^  +  ^  =  0 

dt  duj 


(4) 


where  and  —  ujN^  are  the  density  of  waves  and  the  density  of 
energy  in  the  frequency  space,  correspondingly. 

The  number  of  waves  in  frequency  space  is  connected  with  the  number 
of  weaves  in  wavenumbers  space  as 


/ 


Nuidu  = 


/ 


Tij^dk 


and,  hence,  in  isotropical  case 


Nu,  =  inukk 


^_idk  .  ujc 


,  =  Att 

dw  aC 


-1 

—nk 


a 
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The  energy  flux,  according  to  (  4)  is 


=  —  y  ujN^du:  ==  — 47r 


a‘^C^~ 


UJ 


-3/+1 


ijy) 

-y  + 1 


(5) 


For  y  —  I  (Kolmogorov  spectrum)  the  expression  contains  indeter¬ 
minancy  ^  (see  Sec...)  and  can  be  regularized  using  L’Hopital  rule: 

>1"  dl{y) 


P  =  Att- 


cPC 


dy 


Spectral  density  of  energy  integrated  over  angles  in  /F-space  is  (see  Ap¬ 
pendix  A): 


h  =  dTTCfc  =  A-KR^B{i/Q)k  3 

Comparing  with  definition  of  the  Kolmogorov  constant 


R  —  C kolmP ^  k 


2  , 

3 


we  have 


P'  kolm 


AitR^BIvo) 


pi 


Plugging  in  here  the  expression  (  5)  for  the  energy  flux  P  we  have 


P'kolm 


.1  8  ~22£+3d  r-./  \ 

(47r)3Q:3C3  “  ^^B{vo) 

(T? 


(6) 


For  the  Kolmogorov  spectrum  we  have  y  =  yo  =  1.  Remembering  that 
for  scaling  spectrum  we  have 


y  = -7  + 3a; -  3  =  1 
a  =  ^  —  ax  d 

2/3  A- 3d  , 

7  = - 4 

a 

V  =  xa 
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we  obtain  the  values  of  x,  a  and  p  corresponding  to  Kolmogorov  solution 


(7) 

(8) 
(9) 


It  is  easy  to  see  now  that  power  of  C  for  these  values  of  x  and  q  in  the 
expression  (  6)  equal  0,  i.e.  C  droppes  out  from  the  final  expression  for  the 
Kolmogorov  constant: 


Ckolm  — 


(47r)3Qg.g(i/o) 

(  Mvl 

V  vo=i 


(10) 


where 


dy 


yo=l 


-2^5(1,  ,  ^2"" ,  ^3  “ )  (666)“"°  ( 1  +  i?  -  -  e?" ) 

(<^1  In  1^1  “  ^2  In  <^2  —  ifs  In  <^3)  ^  (1  +  —  ^2  —  ^3)  d(^id^2d^3 

Using  the  expression  (  13)  for  S  ^1 ,  ^1“° ,  ^2''°  5  ^3'*°  ^  from  Appendix  B 


we 


finall}^  get 


m 

dy 


yo=l 

13 


=  (11) 

-# -#)  (61n6-61n6-61n6)x 

6  (1  +  —  ^2  ~  ^3)  <i^id^2d^3 

Expressions  (  10),  (  11)  and  (12)  present  final  formulae  for  a  calculation 
of  Kolmogorov  constant. 

Taking  the  integral  (11)  numerically  we  get 


Ckolm  —  4.7 
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3  Appendix  A 

In  the  Appendix  A  we  are  going  to  calculate  the  spectral  density  of  energy 
Cfc  on  Kolmogorov  solution  (  2),  n  =  no  =  ^. 

The  spectral  component  of  velocity 

Vk  =  J  +  k2)a*^a^Jk-i,dk2 


hence 


1  I/- 

=  2  <  ^k^k'  >=  2  J  ‘^%k2''^kA'^k2^^^  ~  ^'1  +  k2)dkidk2 

where 

<  a^a*-  >=  nj:6{k  -  k') 

Here  brackets  <>  mean  an  averaging  over  an  ensemble. 

It's  easy  to  see  that 


1 


ki  —  k2 ,  ki ,  /f2 


1  1 


kik2 


where 


(2x)5  \k^-k2\^  2(2;r)t 


k  —  ki  —  A* 2 


-  J-  _ 


k^ 


and 


A  .  ^  i_J__ 
^1*2  4  (2%)^ 


,2  1  L‘2\  7,2  (^1  HY 


2{kf  A  k^^)-  P  - 


k^ 


On  the  scaling  solution  (  2)  we  have 


1 


8  {2Try 


I 


nki + ki)  -  e  - 


(^1^2)  ^d(k  —  k\  +  k2)dk\dk2 
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Now  we  will  average  over  the  angles  in  /?-space  This  procedure  is 
reduced  to  the  averaging  of  6-function  over  these  angles.  It  is  easy  to  see 
that 

2  Vi 

<  8{k\  —  ^2  -|-  /js)  >=  ,  ,  , 

/Cl  ^2^3 

where  the  brackets  <>  mean  an  averaging  over  angles  in  A'-space  and  we 
have 


where 


tk  —  R^B{i')k 


5-2i/ 


B{i^)  =  [2(u^  -h  -  1  -  {u^  +  v'^f]  {uvf 


dudv 


Here  the  integration  is  taken  over  half-infinite  strip  A  :  {u  <  u  -f-  1,  u  > 
u  —  l,u  >  u  -h  1}  (see  Fig. 4a).  For  the  numerical  integration  purposes  it’s 
convenient  to  turn  the  system  of  coordinates  by  ^  Coordinate  transformation 


u  = 

V  = 


X  —  y  1 

x^y  \ 

y/2  2 


maps  the  domain  A  into  domain  A'  :  <  y  <  >  0}  with 

Jakobian  =  1 

This  mapping  gives  for  Kolmogorov  spectrum  (//  =  1) 


^C-t)  =  (TvI  +  .^yll 

V3/  {2rY  Ja"  [x -y  +  3  (x  +  y  + 

Here  the  integration  is  carried  over  one  half  of  domain  A',  i.e.  A"  :  {0  < 
X  <  oo,0  <  y  <  (see  Fig.4b). 

From  numerical  integration  we  get 

B  (y )  =  0.047  (12) 
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4  Appendix  B 

In  the  Appendix  B  we  are  calculating  an  analitical  expression  for  the  average 
of  the  matrix  element  over  eight  angles  in  A'-space. 

The  configuration  of  four  vectors  fc,  Ai,  A^,  A3  is  invariant  with  respect 
to  rotations  around  the  vector  p  =  Ai  —  A3  =  A4  —  A2  (see  Fig.l)  as  a  whole. 
This  fact  allows  to  reduce  the  integration  over  four  angles  (j)  to  two  angles  d. 

The  angles  6  can  be  expressed  through  the  absolute  values  of  vectors 
A  =  |A|,  Ai  =  |Ai|,  A2  =  IA2I,  A"3  =  |A'3|,  p  =  |A  —  A2I  =  IA3  —  A2I  as 


cos  6 
cos  0i 
cos  ^2 
cos  03 


A^  +  p^  —  AI 
2pA 

A|  -  p^  -  A| 

2pA’i 

P-P^-H 

2pk2 


H  +  P^-  Pi 


2pA3 


JL  X  X 

Using  these  properties  we  get  the  expression  for  function  ^(l,  ,  ^2“  ?  ^3“ ) 

which  we  will  need  for  our  purposes: 


where  function  R  is  the  following: 


R{ki,K2,K3)  = 

2  ((k2  -  1)^  -  -  (1  +  «2)^)((«3  -  Kif  -  P){P  -  (/Cl  + 

Ji~K2 

{{k3  -  1)^  -  P){P  -  (1  +  «3)^)((/C2  -  /Ci)^  -  At^)(/C^  -  (aCi  +  ^^2^) 

Jl-Kz  K'* 


-  4 


C2P  -C3  +  CiArcTan 


K=K1+«3 


K— 1  ~«2 
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where 


Co  =  l  +  +  + 

r  =  ~  6)^(6  +  ^3)(^2*f3  +  666  +  '^3) 

^  "  U?  +  66  +  ^l)(eH6  +  i) 

C2  =  (1  —  6)^(Ci  +  6^2  +  +  6  +  1) 

C3  =  3(1  -  6)^(1  -  6)^(6 + 6)(^2e3 + ^166 + en 

Q  _  (1  ~  '^2)(1  —  ^3)(^l  +  ^1^2  +  ^2)(^3  +  ^3  +  1)(^2  +  ^2  +  1)(^3  +  6^3  +  ) 

'  ”  \/3(f2  +  isMUl  +  6653  +  6) 

Here  k,  =  ^  and  =  Kf  are  dimensionless  wavenumbers  and  frequencies, 
correspondingly,  2  =  1,3. 
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Abstract 

An  ensemble  of  weakly  interacting  capillary  waves  on  the  surface 
of  ideal  fluid  forms  a  kind  of  wave  turbulence.  Similary  to  the  case  of 
’’classical”  turbulence  in  incompressible  three-dimensional  fluid,  the 
main  physical  process  in  this  turbulence  is  the  energy  transport  in 
A-space  from  large  to  small  scales. 

According  to  the  weak-turbulent  theory,  the  spatial  spectrum  of 
elevations  1^  =<  IVkl  ii^  this  turbulence  is  the  Kolmogorov-type 
spectrum  ^  k~^ .  So  far  this  result  was  not  confirmed  neither 
experimentally  nor  numerically. 

We  developed  a  numerical  scheme  for  direct  simulation  of  the  sur¬ 
face  of  ideal  fluid  based  on  the  expansion  of  the  Hamiltonian  of  the 
surface  up  to  terms  of  fourth  order  and  observed  the  predicted  Kol¬ 
mogorov  spectrum. 


1 


Introduction 


An  ensemble  of  weakly  interacting  waves  in  a  dispersive  medium  can  be 
described  statistically  even  being  very  far  from  the  state  of  thermodynamic 
equilibrium.  Due  to  small  value  of  of  nonlinearity  the  infinite  system  of 
equations  for  corellation  function  in  this  case  can  be  truncated  by  a  consistent 
w-ay  and  reduced  to  one  kinetic  equation  for  ’’wave  numbers”  (or  wave  action) 

3ti 

-^+2^^n^  =  st{n)  (1) 

(see,  for  instanse,  [2]).  Here  7^  is  the  wave  damping  (  or  the  w'ave  pumping 
if  7^  <  0),  st{n)  is  the  ’’collision  term”  corresponding  to  wave  equation. 

The  collision  term  describes  ”cascades”-transport  of  wave  energy  in  K- 
space  to  small  scales  region  (direct  cascade)  and  to  the  large  scale  region 
(inverse  cascade).  The  last  one  existy  only  if  the  total  number  of  waves 
N  =  j  n^dk  is  the  integral  of  motion. 

The  equation 


=  0 

besides  trivial  thermodynamic  solution  has  Kolmogorov-type  solutions 
describing  cascades.  In  a  medium  without  a  characteristic  length  they  are 
powerlike  functions 


nf  ~  k-  (2) 

The  theory  of  weak-turbulent  Kolmogorov  spectra  now  is  far  advanced. 
But  direct  experimental  confirmation  of  these  spectra  are  very  poor.  One  can 
consider  more  or  less  well  confirmed  existence  of  the  Komogorov  spectrum 
for  the  direct  cascade  of  gravitational  wave  on  the  surface  of  incompressible 
deep  fluid 


lui  —  Ct 


(Here  /(u;)-spectral  density  of  surface  elevations,  u-  wave  frequency,  g 
is  gravity  acceleration,  v-wind  velocity,  a  -  dimensionless  constant).  This 
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spectrum  was  theoretically  derived  by  Zakharov  and  Filonenko  [4]  and  ex¬ 
perimentally  observed  by  Toba  [5]. 

Another  way  to  check  the  weak-turbulent  thory  is  numerical  simulation. 
Some  valuable  results  were  obtained  by  numerical  solution  of  kinetic  equation 
(  1)  [7],  [8]-  But  the  kinetic  equation  (  1)  is  itself  a  subject  for  careful 
examination.  Its  derivation  assumes  that  the  phases  of  all  interacting  waves 
are  random  and  are  in  the  state  of  chaotic  motion.  The  validity  of  this 
assumption  is  no  clear  apriory. 

The  right  way  to  check  the  weak-turbulent  theory  and  its  prediction  is 
numerical  simulation  ’’from  the  first  principles”,  i.e.  direct  solution  of  the 
nonlinear  dynamic  equation  governing  propagation  and  interaction  of  the 
waves. 

In  real  cases  these  equations  are  of  two  or  three  spatial  dimensions,  and 
its  numerical  solution  is  not  a  simple  problem.  It  was  done  so  far  for  ‘ID 
Nonlinear  Schroedinger  Equation  [10],  but  in  this  particular  case  Kolmogorov 
spectra  don’t  exist. 

In  this  paper  we  present  results  of  numerical  simulation  of  capillary  waves 
on  the  surface  of  the  incompressible  infinitely  deep  fluid.  In  this  case  only 
direct  cascade  of  energy  takes  place.  Corresponding  Kolmogorov  spectrum 
for  the  surface  elevation  has  a  form  Ik  ~  k~ .  We  will  show  that  this 
theoretical  prediction  is  confirmed  by  direct  numerical  simulation  with  good 
accuracy.  The  developed  numerical  approach  can  be  used  for  solution  of  a 
wide  class  of  problem  pertained  to  interaction  of  surface  waves  and  -  more 
generally  -  other  types  of  waves  in  nonlinear  media. 


Theoretical  background 

We  study  the  potential  flow  of  ideal  incompressible  deep  fluid  with  the  free 
surface.  Let  r  =  {x,y)  is  the  shape  of  the  surface,  ^(r,t)  is  the 

velocity  potential  $  =  $(r,2),  v  =  V$,  A$  =  0,  evaluated  on  the  free 
surface:  It  is  known  [9]  that  under  these  assumption 

the  fluid  is  a  Hamiltonian  system: 

dt  6^  ^  ’ 
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^ 

dt  Srj 


(4) 


Here  H  is  the  total  energy  of  the  fluid  consisting  of  the  kinetic  and  the 
potential  components 


H  —  Hpot  +  Hkin  (5) 

where 

^Po*  =  a|[(l  +  (Vr,f)5-l]dr 
Hkrn  =  \j  dr  L 

Hex'e  a  is  a  coefficient  of  surface  tension. 

Direct  numerical  simulation  of  the  system  (  3),  (  4)  provides  solution  of 
the  boundary  problem  for  the  Laplace  equation  on  each  step  in  time.  In  full 
3D  case  it  is  enormously  hard  problem.  To  solve  the  problem  one  can  use  an 
expansion  in  powers  of  nonlinearity.  For  Fourier  transforms  this  expansion 
up  to  the  quadratic  terms  has  a  form: 


H 

Ho 


H, 


H2 

^k,k2 


^2  ^3  ^4 


Hq  +  Hi  +  H2  +  ... 

~2  x2t  J  +  ^2  +  kz)dhdhdh 

4^2jr)2  j  ^^kklki^y  ki'^k2'^kz^^k  +  +  ^2  + 

k\k2  +  I  1 1  ^2 1 

1^1  II^2|[2(I^1  +  ^sI  +  1^1  +  ^'4!  +  1^2  +  +  1^2  +  ^4!) 


Corresponding  dynamic  equations  are 


^  -  \k\  [[|X-|V’],  X  T]r^^  +  |A:|  \k\  [[|fc|V> 


X  7?,-  ,  X  f], 


dtpf 

~W 


1  ^ 

H — 

2 

O-Ar^r'  + 


\k\ll) 
1 


X  Vf 


.+  [Ar-X  7/? 


—Axp  X  [fcjV) 


-(VV’r4-  [\k\rP 

X  7]  f-\-  Dj^  +  Jp; 


-l^i 


(6) 

(7) 


We  added  to  the  equation  (  7)  a  phenomenological  damping  term  Dr  and 
the  external  force  Fr(t). 

In  the  linear  approximation  the  equations  (  6)  and  (  7)  describe  capillary 
waves  with  the  dispersion  relation 


U>k  =  {(Tk^)2 


One  can  introduce  normal  amplitudes 


According  to  the  weak-turbulent  theory  the  pair  corellation  function 


<  a^a*-,  >=  nj^6{k  —  k') 
satisfies  the  kinetic  equation  (  1),  where 


st{n)  =  j 


^kki  *2  ^ki  ^k2  kki 


^kki  kl 


^kki  k2 


=  47r  I kf6ik-h-k2  -  LOj^^  -  ) 


kki  k2 
^k^^k2  -  ^k'^k. 


^k^k2 


2t2;r 


{u>k,^^k2‘^kD^ 


A 

kik2 


kki 


kko 


In  isotropic  medium  containing  no  characteristic  length  the  dispersion 
relation  is  a  power-like  function 
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UJk  — 


as  far  as  is  a  homogeneous  function 

^tk^eki,tk2 

In  this  case  the  equation 


st{n)  —  0 

has  exact  poAverlike  solution 


rik  = 


c 


pi 

k^d 


(d  is  dimension  of  space),  which  is  Kolmogorov-type  spectrum  describing 
constant  flux  of  energy  in  /\ -space  from  large  to  small  scales.  P  is  the  value 
of  the  energy  flux,  C  is  an  absolute  constant.  For  capillary  waves  a  =  |, 
/?  =  |,  d  =  2.  Hence 


Uk 


C 


pi 

w 


For  the  corellation  function  of  elevation  one  get 


(8) 


h  =<  |7/a-P  >= 


iOk  Ca  Ipl 


This  result  was  obtained  first  by  Zakharov  and  Filonenko  [4].  The  solution 
(  8)  is  linearly  stable  in  framework  of  the  kinetic  equation  (  1)  (see  [11]). 
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Numerical  simulation 


We  realized  numerical  simulation  of  the  system  (  6),  (  7).  In  spite  of  the 
fact  that  the  matrix  element  of  the  kinetic  equation  is  expressed  only 

through  coefficients  of  the  Hamiltonian  Ho,  Hi  we  prefer  to  keep  the  next 
term  H2  in  the  expansion  of  the  Hamiltonian.  The  reason  for  that  is  the 
following  -  it  can  be  shown  that  the  dynamical  system  generated  by  the 
Hamiltonian  Ho  +  Hi  becomes  ill-posed  at  very  low  levels  of  nonlinearity. 
Meanwhile,  including  into  consideration  the  next  term  of  expansion  improves 
the  situation  essentially  (details  of  the  consideration  will  be  published  sepa¬ 
rately).  Moreover,  the  developed  scheme  after  a  minor  modification  can  be 
used  for  numerical  simulation  of  the  gravitational  waves. 

The  equations  (  6)  and  (  7)  are  not  differential  in  the  X-space.  Besides 
taking  derivatives  they  include  taking  of  the  operator  \k\  ((— A)2  in  X-space). 

So,  the  system  can  be  reduced  to  the  set  of  6  PDE  for  the  variables 
interconnected  by  the  consequent  application  of  the  operator  (—A)  2.  It 
allows  to  apply  for  solution  of  the  system  (  6)  and  (7)  the  spectral  code 
using  the  Fast  Fourier  Transform  on  each  step  of  time.  Omitting  the  details 
of  numerical  scheme,  we  reproduce  now  only  the  final  result  of  calculation. 

For  numerical  integration  of  the  equations  (  6)-(  7)  we  used  the  functions 
F  and  D  defined  at  Fourier  space  through  the  following  relations  : 


Fj:  = 

Dk  =  Ik^k 

^ic  ~  ■^(^)) 

=  /  -d^l-l^^ol)'  iik>ko 

[0  0  <  k  <  ko 


Pumping  force  Fk  is  ’’almost”  in  resonance  with  local  linear  frequency  ujk 
of  the  corresponding  Fourier-garmonics,  i.e.  frequency  0^  slightly  fluctuate 
around  exact  value  of  u>k  due  to  small  random  in  time  addition  R{t).  The  form 

/l£H£iy 

of  pumping  amplitude  was  chosen  to  be  axially-symmetric  fk  =  foe^  /  . 

Value  ko  defines  starting  point  of  ’’hyperviscosity”  we  used  in  our  exper¬ 
iments  to  provide  wide  enough  intertial  interval.  Calculations  were  carried 
out  on  the  grid  256  x  256. 
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System  was  driven  by  forcing  F  localized  at  small  wavenumbers.  After  a 
while,  we  observed  formation  of  the  stationary  spectrum  of  waves  carrying 
constant  in  time  energy  flux  to  the  high  k  due  to  nonlinear  interaction  of 
waves.  Observed  spectrum  was  characterized  by  angular  isotropy  and  being 
averaged  over  the  angles  in  A'-space  apeared  to  be  in  a  good  agreement  for  the 
inertial  interval  with  the  spectrum  predicted  by  the  weak-turbulent  theory 
(8). 

This  stationary  spectra  w'ere  obtained  for  and  nonlinearity  levels  ~ 
0.05  (here  Hi  and  are  linear  and  nonlinear  part  of  the  energy). 

Conclusions 

Summarizing  the  results  we  can  conclude  that  the  direct  numerical  simulation 
of  the  dynamic  nonlinear  equation  confirms  an  existence  and  important  role 
of  the  weak-turbulent  Kolmogorov  spectra  at  least  in  the  case  of  capillary 
waves.  Indirectly  this  result  confirms  the  validity  of  the  kinetic  equation  for 
a  description  of  the  wave  turbulence.  We  hope  that  the  developed  effective 
approach  will  allow  us  to  study  numerically  other  types  of  wave  turbulence, 
first  of  all  -  the  behavior  of  a  system  of  wind-driven  gravitational  waves  on 
the  sea-surface. 
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